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Abstract —Future wireless communication systems are de¬ 
manding a more flexible physical layer. GFDM is a block 
filtered multicarrier modulation scheme proposed to add multiple 
degrees of freedom and cover other waveforms in a single 
framework. In this paper, GFDM modulation and demodulation 
will be presented as a frequency-domain circular convolution, 
allowing for a reduction of the implementation complexity when 
MF, ZF and MMSE filters are employed in the inner and 
outer receiver operation. Also, precoding is Introduced to further 
increase GFDM flexibility, addressing a wider set of applications. 

Index Terms —precoding, low-complexity, multicarrier modu¬ 
lation, GFDM. 

I. Introduction 

Fifth generation (5G) networks are requiring a new level 
of flexibility on the physical layer (PHY). Although higher 
throughput keeps pushing the spectral efficiency to higher 
standards, new services will also demand very low latency, 
massive capacity for multiple connections and very low power 
consumption. The PHY must be flexible to address several 
different scenarios. 

Generalized Frequency Division Multiplexing (GFDM) ||T| 
is a recent waveform that can be engineered to address 
various use cases. GFDM arranges the data symbols in a 
time-frequency grid, consisting of M subsymbols and K 
subcarriers, and applies a circular prototype filter for each 
subcarrier. GFDM can be easily configured to cover other 
waveforms, such as Orthogonal Frequency Division Multiplex¬ 
ing (OFDM) and Single Carrier Frequency Domain Multiple 
Access (SC-FDMA) as comer cases. The subcarrier filtering 
can reduce the out-of-band (OOB) emissions, control peak 
to average power ratio (PAPR) and allow dynamic spectrum 
allocation. GFDM, with its block-based structure, can reuse 
several solutions developed for OFDM, for instance, the con¬ 
cept of a cyclic prefix (CP) to avoid inter-frame interference. 
Hence, frequency-domain equalization can be efficiently em¬ 
ployed to combat the effects of multipath channels prior to the 
demodulation process. With these features, GFDM can address 
the requirements of 5G networks. 

The main disadvantage of a flexible waveform is the com¬ 
plexity required for its implementation. Reducing the GFDM 
complexity is essential to bring its flexibility to 5G networks. 
The understanding of GFDM as a filtered multicarrier scheme 
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leads to a modem implementation based on time-domain 
circular convolution for each subcarrier, performed in the 
frequency-domain |[T]. This paper follows the well known 
polyphase implementation of filter bank modulation Q-Q 
and applies it to GFDM. 

As one main contribution, the investigation in this paper 
reveals that the Poisson summation formula © can be utilized 
in the demodulation process, considering a per-subsymbol 
circular convolution and decimation in the frequency-domain. 
This operation can be performed as an element-wise multi¬ 
plication with subsequent simple M-fold accumulation in the 
time-domain. This simple re-orientation of the data symbol 
processing allows for a considerable reduction in complexity, 
which makes GFDM applicable in a wide range of scenarios. 

This paper also contributes to reformulate the modula¬ 
tion and demodulation as matrix-vector operations based on 
M{K)-pomt discrete Fourier transforms (DFTs). As an out¬ 
come, it becomes evident that DFT can actually be considered 
as a precoding operation. Moreover, they can be replaced 
by other transformations as well. This observation facilitates 
arbitrary precoding of the GFDM data, which adds a new level 
of flexibility to the system. 

It is worth mentioning that the proposed low complexity 
signal processing complements the work in © because the re¬ 
quirement for block alignment is loosened. It can be usable for 
supporting pipeline inner receiver implementation, particularly 
when building synchronization and channel estimation circuits 
for embedded training sequences. It can also be beneficial to 
develop future non-linear and recursive detection algorithms. 

H. Classical GFDM Description and 
Low-Complexity Reformulation 

Consider a wireless communication system that transmits 
data in a block-based structure that consists of K subcarriers 
and M subsymbols. Let N = KM be the total number of 
data symbols in the block and dk,m, k = 0,..., K — 1, 
m = 0,..., M — 1, denote the complex valued data symbol 
that is transmitted on the fcth subcarrier and mth subsymbol. 
The classical description of GFDM signal generation © is 
given by 

K-l 

= X ®M'n] (1) 

M-1 

with dk[n] = ^ dk,m5[n - mK], (2) 

m—0 
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where n = 0,..., — 1, ® describes circular convolution 

carried out with period N and g[n] denotes the impulse 
response of the transmit prototype filter. 

Using the A^-point DFT Q can be carried out as 

K-l 

x[n] = '^{FN{g[n]eyi.-p{j2T:^n)} FN{dk[n]} }, (3) 

where (•)^ denotes the remainder modulo N and 
G[f] = J-N{g[n]}. Due to periodicity of Tm with period M 
and / = 0,... ,MK—1, J-M{dk,( )} concatenated K times. 
This leads to the low-complexity modulator implementation 
described in 0 - Now, expressing the convolution in 0 
explicitly, the transmission equation can be reformulated as 

0 

M-1 K-l 

x[n]= ^ g[{n - mK )^ dk,m exp{j2n-^n). (4) 

m—0 k—0 

' V-^ 

The sequence obtained from the Inverse Discrete Fourier 
Transform (IDFT) is concatenated M times in time-domain 
due to n = 0,..., MK — 1. Such multiplication in the time- 
domain expresses the convolution of the subcarrier filter with 
the data in the frequency-domain for each subsymbol. 

Assuming a flat and noiseless channel, the convolution of 
the received signal y\n] = x[n] with a demodulation filter 
7 [n], is given by 

dk,ra = ( 7 * [(-«) at] ® 2/W exp(-j27r;^n)) \n=7nK- (5) 

With the help of DFT, 0 is written as 
dk,m = exp(-j27r An)}}|™x, 

'-V--V-^ 

r*[/l y[</+fcM)jvi 

(6) 


where J^Ar{ 7 [n]} = r[/] and J-N{y[nW = y[f]- For example, 
this leads to the low-complexity demodulator implementation 
described in 0 for the Matched filter (MF) 7 [n] = g[n\, which 
was later extended to more general filter types in Q. Now, 
directly expressing the convolution and sampling in^^ leads 
to the time-domain multiplication 

N-l 

dk,m = l*[{n - mK)^]y[n] exp{-j2n^n) 

= - mK)j^]y[n]}\kM, 

where the A^-point DFT is only evaluated at every Mth 
sample. According to the Poisson summation formula 0, this 
can be reduced to a ^ = K point DFT 

{ M-l 

m=0 

Eq. 0 describes a sampled short-time Fourier transform 
(STFT). This is clear since the GFDM demodulator is actually 
performing a Gabor transform 0, which can be understood as 
a sampled STFT. For this reason, 0 does not allow to employ 




Fig. 1: Overview of block structure and terminology. 


different prototype filters per subcarrier. This constraint is 
circumvented by performing equalization before demodulation 
or the presented scheme can be used to initially access the 
channel state information by demodulation of training data. 
Also, the proposed scheme facilitates the implementation of 
odd number of subsymbols 0 avoiding the use of M-point 
DFT and IDFT required in |[^ 

Notice that 0 holds when MF, zero-forcing (ZF) or 
minimum mean square error (MMSE) filters are used. ZE 
or MMSE filters can deliver good symbol error rate (SER) 
performance, while ME needs successive interference cancel¬ 
lation (SIC) to remove the SER floor caused by the self¬ 
interference. Hence, 0 shall be used to implement ZE and 
MMSE demodulators, while ME can be considered when the 
transmit Alters are orthogonal or Offset Quadrature Amplitude 
Modulation (OQAM) is employed. 

III. GEDM Matrix Model 
Arranging the data symbols dk^m in ^ two-dimensional 
structure leads to the data matrix D = (dg di • • • dM-i) , 
where d^ = (rfo.m f^i.m • • ■ denotes the data 

symbols transmitted in the mth subsymbol. 

According to the GEDM description in Section the 
rows and columns correspond to the time and frequency 
resources in Eig. 0 respectively. Hence, the A:th row of D 
represents the data symbols transmitted in the A:th subcarrier, 
d/c = {dk,o dk,i ■ ■ ■ dk,M-i) ■ The time-domain circu¬ 
lar convolution of the modulation process can be expressed in 
the frequency-domain as 0 

^ Y (9) 

k=0 

Eor each subcarrier, the data vector is taken to the 
frequency-domain by the M-point DET matrix The 

corresponding time-domain up sampling operation is real¬ 
ized in the frequency-domain by duplicating the transformed 
data symbols vectors K times, using the repetition matrix 
where Im is an M size identity matrix, 
j is a 2 X j matrix of ones and ® is the Kronecker product. 
Each subcarrier is then Altered with G = diag(WArg), where 
diag(*) returns a matrix that contains the argument vector on 
its diagonal and zeros otherwise and g is the vector containing 
the transmit Alter impulse response. An up-conversion of 
the fcth subcarrier to its respective subcarrier frequency is 
performed by the shift matrix 

p(G = ^ ^ (10) 
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Fig. 2: Block diagram of low-complexity modulator and demodulator with precoding. 


where returns the circulant matrix based on the input 

vector and is a column vector where the fcth element is 
1 and all others are zero. The K subcarriers are summed and 
transformed back to the time-domain with to compose 
the GFDM signal. On the demodulator side, the recovered data 
symbols for the A:th subcarrier are given by 

dfc = ^WH ^ r ^ W^y, (11) 

where y is the equalized vector at the input of the demodulator, 
r = diag(W N'y) with 7 being the demodulation filter impulse 
response, e.g., MF, ZF or MMSE filters. 

The modulation and demodulation processes can be simpli¬ 
fied just by changing the processing order of the data symbols, 
as derived in Q. In this new matrix model, the GFDM vector 
is given by 


M-l 




( 12 ) 


m—0 


Observe that (12i is similar to a polyphase filter structure, but 
with the difference that cyclic time-shifts are used here and 
frequency-domain convolution is performed in time-domain 
as element-wise vector multiplication. With this approach, 
the first step is to obtain a time-domain version of by 
multiplying it with an inverse DFT (IDFT) matrix W^. M 
times upsampling in frequency-domain is performed in time 
by duplicating the transformed data symbols with a repetition 
matrix Each subsymbol is then pulse-shaped with 

g. p(”^) shifts the mth subsymbol to its respective position in 
time. The GEDM signal is obtained by summing all the pulse¬ 
shaped subsymbols, with no need of the A^-IDET domain 
conversion in @- 

The demodulation operations are also simplified by using 
the circular-convolution in frequency-domain, leading to 


= Wk 


diag( 7 )(p('"))^y. (13) 


While the new approach describing GEDM modulation 
and demodulation as convolution in the frequency-domain 
reduces the implementation complexity, in the next section 
this principle will be expanded to a more general structure, 
where domain conversions will be understood as a precoding 
process. 


IV. Expanding GEDM features with precoding 

Eqs. (12 1 and ( [T3] l bring a new interpretation of the GEDM 
chain that further increases the flexibility of this waveform. 


Let the transmit vector from 0 be redefined 


M-l 


X = ^ p(m)diag(g)R(M,if)^^ 


m—0 


where 


Srr,. = Td„ 


(14) 

(15) 


is the information coefficient vector to be transmitted and 
T is a transformation matrix. Eor classical GEDM, where 
a time-frequency grid is used to transmit the information, 
T = Wj^. However, (15 1 can be seen as precoding the 
information vector by a generic matrix T. Different matrices 
can be used to achieve different requirements. Hence, by 
transmitting precoded data, GEDM can be even more flexible 
to address requirements of future wireless networks. 

On the receiver side, the recovered coefficients are given by 
T / , T 


Sm = 


j'R(M.ir)^ diag(7)(p(" 


y- 


(16) 


The estimated data symbols are obtained by reverting the 
precoding operation on ( [T6| l, leading to 

dm = SSm- (17) 


where S removes precoding. Eor conventional GEDM, 
S = Wif. A block diagram that described the proposed 
GEDM modem is presented in Eig. The demodulation part 
can be implemented to run continuously, which can be used 
for tracking specific properties of the estimated information 
coefficient vector Sm containing training data. 

The precoding concept can be enlarged to embrace any 
transform taking into account the two dimensions of the 
GEDM data block D, which is a bi-dimensional structure 
defined as a frequency-time grid spanning over several sub¬ 
carriers and subsymbols, as depicted in Eig. [T] Hence, in 
the default configuration, the columns of the matrix contain 
data defined in the frequency-domain, while the rows have 
data in the time-domain. The precoded data block A can be 
generalized to 


A = Te • D • Tr , (18) 

KxM KxK KxM MxM 

where D is now a data block defined in an arbitrary domain, 
Tc is a AT by iT precoding matrix applied to the columns and 
Tr is a M by M precoding matrix applied to the rows of 
D. Eor instance, choosing Tc and Tr to be Eourier matrices 
yields four domains for the GEDM data block, as given in 
Tab. [I] , where the domain frequency-time (ET) corresponds to 
standard GEDM. The choice of the data domain has impact 
on the robustness of the system regarding time and frequency 
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Fig. 3: OOB, PAPR and complexity analysis for precoded GFDM. 

FT/FF; M = 16, Mo„ = 15, AT = 128, ATon = 75, TT/TF: AT = 1, M = 1200, OFDM; N = 2048, = 1200. 


TABLE I: Precoding in time (T) and frequency (F) domains 


Domain 

Operation 

Number of multiplications 

ET 

A 

= WH Dlif 

MK\og^{K) 

TT 

A 

= IkDIm 

0 

EF 

TF 

A 

A 

= WH DWH 

= IkDWH 

MK log2 {K) +KM log2 (M) 
KM log2(M) 


selective fading. The concept is not limited to DFT and, in 
general, any meaningful transform can be applied to the rows 
and columns. 

Taking the concept one step further, ( [T^ can be formulated 
as 

vec(A) = vec(TcDT,). (19) 

With the help of the Kronecker product, the right part of ( [T9| ) 
can be expressed as 

vec(T,DTr) = (t;;^ 0 T,) vec(D). (20) 

The resulting matrix 

r=(TTGTe) (21) 

can be seen as an even more general way of precoding that 
is not restricted to subcarriers or subsymbols but allows 
arbitrary coupling between any elements of d = vec(D). The 
precoding is then applied leading to ^ = Td. 

V. Out of Band, PAPR and Complexity Analysis 

In particular, the FT mode allows for use of null sub¬ 
symbols and subcarriers to support non-continuous bandwidth 
application, achieving a low OOB radiation as shown in 
Fig. Although this is not possible with the time-time (TT) 
case, where the system essentially becomes a single-carrier 
system occupying the whole bandwidth, this precoding is 
especially beneficial for power-limited systems that cannot 
afford complex precoding operations with a low PAPR, as 
the PAPR can be significantly reduced. In contrast, frequency- 
frequency precoding (FF) results in a DFT-precoded GFDM 
system, with a spectrum similar to OFDM. In FF mode, data 
is defined directly in the frequency domain, according to ([^, 
with one particular subsymbol only located on at most two 
frequency bins, assuming a transmit filter that is bandlimited to 


two subcarriers. This property simplifies equalization, as only 
ICI remains, but on the other hand reduces frequency diversity 
in channels with strong frequency-selectivity. The last (TF) 
case also corresponds to the single carrier, but with symbols 
spread in time, creating potential for exploiting time-diversity. 

The complexity of the precoding operation of the GFDM 
transceiver based on DFT algorithms are evaluated in Table 
Number of complex valued multiplications, a costly operation 
in implementations, is the chosen criterion for the analysis. 
The comparison is carried out under the assumption that 
N = KM complex data symbols are transmitted. As the base¬ 
line for the comparison in Figj^ MN multiplications from 
( [Tb] ) are added to the precoding operations and the OFDM 
algorithm requires naturally the complexity of a DFT opera¬ 
tion, with A^log 2 (W) multiplications. The GFDM modulator 
and demodulator proposed in ||T) require the complexity of 
the M-point DFT algorithm with additional AT times complex 
multiplications over repeated chunks of M complex samples, 
resulting from a DFT operation. Even when the frequency 
response of the demodulation filter is assumed to be sparse, 
e.g. when roll-off is small and a ZE filter span L can be smaller 
than the total number of subcarriers, the modem scheme 
presented in this paper is less complex than the solution 
proposed in ||Tj. 
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